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Using canonical transformations we diagonalize approximately the Hamiltonian of a gaussian wire with 
Rashba spin-orbit interaction. This proceedure allows us to obtain the energy dispersion relations and the 
wavefunctions with good accuracy, even in systems with relatively strong Rashba coupling. With these 
eigenstates one can calculate the non-equilibrium spin densities induced by applying bias voltages across 
the sample. We focus on the ^-component of the spin density, which is related to the spin Hall effect. 
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The spin-orbit interaction opens up the possibility to manipulate the electron spin using electrical means, 
either with applied bias or gate voltage Qj. O ne manifestation of such electrical spin control is the spin 
Hall effect lEJESl where spin currents are created via an interplay between spin-orbit interaction and the 
applied electric field. When the spin Hall effect is considered in systems with finite transverse size the 
edges plays an important role. In disordered systems one would expect a building up of spin density (spin 
accumulation) with opposite sign at the two edges. This picture is not necessarily valid in ballistic systems 
and the wavefunctions (or wavepackets) themselves become important quantities |5|. The wavefunctions 
can then be used to calculate the bias voltage induced spin density 171 [8l l9l ITOl ITT1 IT2l IT3l IT4l IT5l . 

In this paper we calculate analytically an approximate eigenspectrum of the Rashba Hamiltonian in a 
gaussian quantum wire. This allows us to calculate analytically the relevant matrix elements and density 
of states for each transverse mode in the wire with good accuracy, even for relatively strong spin-orbit 
coupling. The Hamiltonian for the quantum wire made in a 2DEG with Rashba spin-orbit interaction, for 
a given wavevector k along the wire, is 

h 2 k 2 1 2 2 a. . 
H = — \--mu)y + —{nka y - py<T x ) (1) 

iLillh — it 

h 2 k 2 + 1 a t a, + . 

+ huj(a'a + — ) + T z ak ^(a 1 r + + ar_) H T - + aT +)- (2) 



2m v 2' ^/2t V2£ 

In Eq. d2J we introduced for convenience new spin matrices where t z = a y and t± raises and lowers 
between the a y eigenstates. For positive momenta only states \k, n, 1) and \k, n + 1, |) cross in energy. All 
other states are separated by at least fku, which allows us to treat the term proportional to (aV + + ar_) 
perturbatively in k so £ = a/lhoj. Note that our Hamiltonian is time reversal symmetric (there is no external 
source of magnetic field). Hence we can focus on the posive k states and use the Krames relations e„. s (k) = 
e n .- s (— k) and ipk,n,s{ r ) — 4 } -k,n,-s{ r ) to obtain the negative k states. 

As mentioned above, for the positive momentum states the perturbative term is 

V = -2^ (T+at+T - a) ' (3) 
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and the non-perturbed Hamiltonian is 

h 2 k 2 / t A i a , t\ 

H = — h ftw(a T a + -) + ar z k + —^(t+cl + r_a T ). (4) 

2m 2 V2£ 

We define an effective Hamiltonian iJ e ff = e s He~ s which results in 

#eff = Ho + i[5,y] + |[5,[5,y]] + ((^) 4 ), (5) 

with S chosen to fulfill the usual condition V + [S, Hq] = 0. The S which satisfies the above requirement 

is 

The new terms generated by 5 in the effective Hamiltonian are 

1 _ a 2 (r,(2 a tq+l)-l) Q 3 (r + aat a + r +a t 3 +h.c.) 

2^ J - Ai 2 (^ + 2ak) + V2AP^ + 2ak) +0 ^ J >>- (7) 

The term proportional to r + at 3 and its hermitian conjugate are perturbative, i.e. they only couple terms 
which are separated by at least 3hw. Furthermore, one can show that the term h[S, [S, V]] in Eq. © con- 
tains no 'diagonal' terms^ and thus the corrections to the Hamiltonian, and consequently its eigenenergies, 
are at most o((fc so £) 4 ). The resulting effective Hamiltonian is thus 

h 2 k 2 _ , + , , , a 2 (r z (2at a + 1) - 1) 
H eS = — — + Kw{a)a + 1/2 + akr z + V ; 
2m 4£(fiw + 2afc) 

a ( ( a 2 a^a \ ( a 2 a^a \ A 

t+0 1 - T7T , „ „ + 1 - — - — — r_a T . (8) 



2V2£ V V 4£hw(nw + 2ak) J \ Athw{hw + 2ak) 
The above _ff e ff can be exactly diagonalized using a generalized rotation matrix in spin space 

cos(G[n + l]/2) sin(( ^±l 1/2) a 

„ Qt sin(eg+l]/2) cos( @ N/2) 



U — I +sinfe[n+ll/2l /^tin ,„n I ' ( 9 ) 



where n = a T a and 

/»2fr (fc so ^) 2 (ft+l) \ 

cor- , m V 2 so 2(i+2fe so «) J 

,// r l /, nil, (fcso^) 2 (ft+l) \ 2 , (fc so l) 2 (n+l) (k w i)Hh+l) \ 2 
Y V 2 2(l+2fc so ^fc) ^ 2 ^ 4(l+2fc so «) J 

The resulting Hamiltonian iJdiag = U'H e gU is diagonal in the spin and ladder operators and its eigenen- 
ergies for positive k are 

^ ^(^4-^*- 2(1^) ) (U) 
^(Jb) = ^(^ 2 +-- 2(1 f^2 fc) +A„( fc )), »>0 (12) 

£n, T (fc) = ^ Ql 2 fc 2 + (n + 1) - y^fc) ~ A »+i(fe)) ■ ™>° (13) 

a /7 1 " {kJ) 2 ^ V , (fc so ^) 2 A (fc so ^) 2 n V „^ 

A " (fc) = n2~ hJk - 2 { l + 2k so i 2 k) ) + — n { 1 - 4(l + 2k s J 2 k) ) ' (U) 



' Here, diagonal refers to a^a, t 2 and also the terms containing r+ a and r_a^ 
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Fig. 1 Energy dispersions of Eqs. dl lb — <|13b are plotted (solid lines) for the first 6 oscillator levels for k so £ — 0.2 in 
(a) and k so l = 0.4 in (b). The dashed lines are the results of numerical calculations. Note that the numerical curves 
fall ontop of the analytical curves except for the highest energy curves in (b). 



The negative k branch of the spectrum is obtained from the Kramers relations. In Fig. Q]the energy dis- 
persion for the full spectrum is plotted, along with results of numerical calculations (dashed curves). Even 
for a relatively large (k so £ = 0.4) perturbation parameter the numerical and analytical curves are barely 
distinguishable except for the highest transverse states, see Fig.QJ. 

Spin densities: Knowing the wavefunctions and eigenenergies allows one to calculate many physical 
quantities analytically, even non-equilibrium ones. As an example we consider spin densities (transverse 
to the wire direction) and focus on the the average value of a z 8{y). For an applied bias voltage, (a z S(y)} 
is related to the spin-Hall effect [16|. Since the Hamiltonian is time reversal symmetric there will be no 
equilibrium spin density in the system [ 17 1. Applying a bias voltage lifts this symmetry since Kramers pairs 
become unequally occupied. For spin densities (and any quantity which is odd under time reversal) only 
states in the bias window need to be considered. Starting from the non-equilibrium Green's functions lfl8l 
it is possible to show that for kT = the spin density is 



( T v S (v))\eV = ^2(ipkns\T v 5(y)\tfj kns ) f (e ns (k) - (fi + eV)) - f(e ns (k) - ^i) 



(15) 



kns 



where the sum is over all k whose velocity is positive (assuming that eV > 0), i.e. the sum contains not 
only k > states but can also include negative momentum states. Although the bias induced spin densities 
are know analytically through the wavefunctions, the equations are quite long and complicated. As an 
example we give the following matrix element 

{^kn]\T x 8{y)\ipkn]) = -sin(6[n+ l])(j) n (y)<p n+1 (y) 

+ i + 2fcl^fc { cos(e[ " + ^My^n-iiy) 

- sin(e [n + l])/2) 2 0„ +1 (y)4> n+2 (y) } , (16) 

where (y\n) = 4> n (y) are the harmonic oscillator eigenstates. The matrix element in Eq. ( fl~6] l has the 
symmetry property (Vwk^^l^fcns) = -{^kns\o z 8{-y)\^ kns ). This is consistent with the spin-Hall 
effect phenomenology where the transverse spin current leads to opposite spin polarization at the two edges 
lfl6l l5l. This property is a direct result of the parity of the harmonic oscillator states 4> n (y)- 

The spin densities (normalized to the wire length) for different values of the chemical potential /i are 
plotted in Fig. [2] assuming linear response eV -€1 fuo. As expected the spin density is odd in y and the 
highest transverse mode contributes the most due to its density of states being the highest. Note that the 
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size of the spin density is comparable to the induced particle density since the matrix element in Eq. (15[ 
is order unity due to the strong coupling of adjacent oscillator states with opposite spin (see Eq. ([8])) 

To summarize, we have calculated the eigenspectrum of a gaussian quantum wire with Rashba inter- 
action. The spectrum allows us to calculate the non-equilibrium spin density analytically. We presented 
results for the a z density for a few low lying oscillator states but the formalism can well handle higher 
states also. 
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